Mechanical limits of viral capsids 
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We study the elastic properties and mechanical stability of viral 
capsids under external force-loading with computer simulations. 
Our approach allows the implementation of specific geometries 
corresponding to specific phages such as 029 and CCMV. We 
demonstrate how in a combined numerical and experimental ap- 
proach the elastic parameters can be determined with high pre- 
cision. The experimentally observed bimodality of elastic spring 
constants is shown to be of geometrical origin, namely the pres- 
ence of pentavalent units in the viral shell. A criterion for capsid 
breakage is defined, which explains well the experimentally ob- 
served rupture. From our numerics we find for the dependence 
of the rupture force on the Foppl-von Karman (FvK) number 7 a 
crossover from 7 2 / 3 to 7 1 / 2 . For filled capsids high internal pres- 
sures lead to a stronger destabilization of viruses with a buckled 
ground state than unbuckled ones. Finally, we show how our nu- 
merically calculated energy maps can be used to extract informa- 
tion about the strength of protein-protein interactions from rupture 
experiments. 

thin shells | membranes | biomaterials 

Bacteriophage capsids have astonishing elastic properties. They 
withstand extreme internal pressures exerted by their densely 
packed DNA. DNA packaging experiments on phage 029 [1| and 
theoretical arguments [3l HI \5\ [6) imply that phage capsids store 
their DNA under ^50atm. This pressure is necessary to inject the 
DNA into the prokaryotic host cell Q. It has been shown experimen- 
tally that the ejection pressure from A-phages is ^20atm [8|. The 
remarkable robustness under high pressures motivated recent nanoin- 
dentation studies [9], which showed that capsids also resist external 
point forces up to ~lnN. 

In contrast, viruses that infect eukaryotic cells, e.g. CCMV and 
polyoma, penetrate their host cell. Their DNA is released by subse- 
quent disassembly of the shell. Correspondingly, these viruses do not 
store their DNA under high pressure. Nevertheless, their resistance 
to external forces is remarkably strong |[T0l . 

Viral capsids are composed of a small number of different pro- 
teins, which cluster to morphological units ("capsomers"). These 
units are put together in a highly regular fashion described by the 
"quasi-equivalence" principle by Caspar and Klug [ 11 1. Due to the 
presence of capsomers the surface of capsids has a discrete structure. 
Mathematically, these capsomers correspond to the vertices of a reg- 
ular triangulation of a sphere. For such a triangulation the number 
of vertices V, connecting edges E, and associated faces F have to 
fulfill Euler's theorem F — E + V — 2. This theorem implies that 
capsids of icosahedral symmetry have 12 pentavalent morphological 
units ("pentamers") embedded in an environment of hexavalent units 
("hexamers"). 

Viral capsids can undergo elastic and bending deformations. Gen- 
erally, the elastic properties of such thin shells (of typical length scale 
R) depend only on a single dimensionless parameter, the FvK num- 
ber 7 = R 2 hie /Kb set by the ratio between elastic modulus K e and 
bending rigidity Kb. It has been shown in Ref. [12], that at a critical 
7b = 154 spherical shells undergo a buckling transition, in which 
the shell acquires a more faceted shape. This transition is caused by 
large strains associated with the pentamers which are reduced upon 



buckling into a conical shape. 

Typical 7 values of viral shells lie in the range from below the 
buckling threshold (e.g. alfalfa mosaic virus with 7 ~ 60 [13']) up to 
several thousand (e.g. phage T4 with 7 c± 3300 Q3)). An extreme 
example is the giant mimivirus (diameter of ~ 600nm) with a Fvk 
number 7 ~ 20000 [15']. Even higher 7-values are possible for arti- 
ficial capsules such as vesicles with crystallized lipid membrane fT6l 
(7 10 5 ) or polyelectrolyte capsules [17] (7 ~ 10 6 ). 

The elastic properties of capsids can be probed in scanning-force 
microscopy (SFM) experiments [9 10]. Depending on the strength of 
loading two regimes are explored. Small forces lead to a completely 
reversible deformation of capsids. Here, the linear and nonlinear 
regime of thin spherical shell elasticity can be studied. In this regime 
the deformations explore the global elastic properties. Larger forces 
(> InN) cause irreversible changes in the shell structure commonly at- 
tributed to bond rupture. Rupturing studies therefore give information 
about the molecular interactions between capsomers thus elucidating 
local features of shell mechanics. 

Viral shells have also been the subject of numerical investigations. 
The dependence of virus shape on the FvK number was analyzed in 
Ref. JT2l . Only recently, the elasticity of capsids has come into the 
focus of numerical studies [18] which are based on a discretization 
scheme introduced for crystalline membranes in Ref. JT9l . In other 
approaches, the shells are directly built up of proteins l20l or cap- 
somers [21 ] with specified (spatially varying) interactions. On this 
basis the stability of capsids against internal pressure has been stud- 
ied (22). 

Here, we generalize a numerical approach developed for the in- 
vestigation of vesicles [23] to viral capsids. The big advantage of our 
method is its applicability to arbitrary geometries. In particular, our 
discretization of the bending energy does not depend on the underly- 
ing triangulation of the viral surface. Therefore, our method produces 
highly stable and reliable results even under high local strain which 
allows us to investigate the rupture of mechanical shells. In our sim- 
ulations, we can systematically vary the elastic moduli and geometry 
of the capsids and probe their mechanical response to external distur- 
bances. We can even implement specific geometries, corresponding 
to specific phages and viruses, and determine (by a direct compari- 
son with experimental indentation experiments) with high precision 
their elastic parameters, such as linear and nonlinear spring constants, 
and the FvK number. In our simulations, we are also able to access 
features which are not observable in experiments. For example, by 
measuring the local strain we are able to determine numerically the 
spatial distribution of rupture probabilities across the capsid surface. 
We will show that experimental deviations from this distribution may 
be used to draw conclusions about the spatial variation of protein 
binding strength. 
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Fig. 1 . The viral structure is represented by a triangulated sphere. The pen- 
tavalent vertices represent the centers of the pentameric units. Vertices are con- 
nected by elastic bonds with spring constant k e . The local mean curvature is 
calculated from the normals of neighboring facets and their shared edge. 

This paper is organized as follows: after a short summary of the 
numerical methods, we first study the global response of capsids to 
externally applied point forces. The numerical results for (empty) 029 
and CCMV are compared with experimental data. Next, we analyze 
the local, irreversible response to indentation forces and we discuss the 
dependence of local rupture-probability on the geometry of the cap- 
sid. Finally, our analysis is extended to filled capsids whose elastic 
parameters and mechanical stability is studied in the last 2 sections. 

Methods 

We have performed numerical minimization simulations of tri- 
angulated surfaces representing the surface of capsids with elastic 
(stretching) and bending energy. In our (small mesh-size) triangu- 
lation every capsomer is represented by several vertices. Therefore, 
these units also have some flexibility and elasticity. Such a discretiza- 
tion is suitable to determine the shape of viral shells 1 12] and yields 
results which do not depend on the number of vertices (in contrast 
to coarse triangulations in which each capsomer is represented by a 
single vertex ['18']). The triangulation represents the underlying icosa- 
hedral (quasi-) symmetry of the capsid. In particular, the pentavalent 
vertices correspond to the centers of the pentamers. 

In our model, the vertices are connected by harmonic springs. 
Any deviation from the preferred inter- vertex distance / gives rise to 
an elastic energy (T9), 

(iJ) 

where vertex i is at position n and the sum extends over nearest 
neighbor pairs only. k e is related to the 2D Young modulus K e via 

The bending energy is given by 

Eb = Kb'%2j* r , [2] 

i 

where the sum extends over all vertices i. Here, Ai is the area assigned 
to the i-th vertex given by 1/3 of the area of all adjacent triangles. Hi 
is the mean curvature of vertex i (with coordination number Zi) given 
by ED 

Hi = 1 gjj • ( n 3,l X n 3,r) [3] 

2 2-^ \n 3 \r\\n 3 \i\ + rij tr • n jt i ' 

with eij = Vi — rj, and rijj and n 3 , r are the non-normalized surface 
vectors of the facets left and right of the edge connecting vertex i and 
j, see Fig. Q] 



Starting from a triangulated sphere the shape of minimal energy is 
found by minimizing the total energy E — E e + Eb using a conjugate 
gradient algorithm l24l . In order to simulate the SFM-experiments 
mentioned in the introduction vertices i at the point of loading were 
moved away from their equilibrium position n. We thus work in 
an ensemble of prescribed indentation X rather than in an ensemble 
of applied force. All vertices are constrained to lie above a (virtual) 
plane at z = 0. The shape of the triangulated surface under these 
constraints is again found by minimizing the total energy. 
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Fig. 2. Force-distance relations of locally deformed viral shells. For small 7 
the force-distance relation is linear F ~ £ for the full range of applied forces (A). 
Above the buckling threshold the structural inhomogeneity of the capsid material 
is reflected by the increasing differences in the elastic response of pentamers and 
hexamers (B) and (C). Above 7 ~ 1000 the inversion transition of pentamers 
(see Fig. 0J3) causes softening, which leads to a decrease in F(£). The influence 
on hexamers is much weaker. At higher 7 pentamers may even snap into a new 
stable inverted configuration, which leads to discontinuities of The dotted 

blue line in (C) represents metastable conformations which have a higher energy 
than stronger deformed ones. Hexamers do not show a discontinuous inversion 
transition. 
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The shapes of bacteriophage 029 and the CCMV plant virus 
have been quantitatively characterized by Cryo-EM and X-ray studies 
l25l l26 l. Phage 029 has an average equatorial radius R ~ 21nm and 
a shell thickness of h ~ 1.6nm. Its FvK number is 7 ~ ll(R/h) 2 ~ 
2000 implying that the shape of 029 is noticeably buckled. For CCMV 
the corresponding values are R ~ 12nm and h c± 2.5nm. Thus, 
7 ~ 300 and CCMV is slightly buckled. 

Results and Discussion 

Spring Constants. As mentioned in the introduction loading with 
small forces probes the global elastic behavior of viral shells. Here, 
we analyze the influence of internal structure on the elastic response 
of capsids by comparing the numerical simulations with analytical 
results for small deformations. Furthermore, we extract the values 
of the elastic moduli of 029 and CCMV by a direct comparison of 
numerical and experimental force-distance relations. 

In a first step we have numerically calculated force-distance 
curves for viruses indented by a distance X = for < £ < 0.4 
(varying with step-size A£ = 1/250). By minimizing the total en- 
ergy the (dimensionless) force F is obtained by F = —AE/ («&A£), 
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Fig. 3. Linear spring constants as a function of 7 for empty (solid lines) and 
filled (dashed lines) icosahedral capsids. The black dashed line represents the 
analytical solution (see suppl. mater.). For 7 < 75 the spring constants of 
the empty capsid follow closely the theoretical prediction. Above 75, the spring 
constant of hexamers ~ 7 1 / 3 . Pentamers become stiffer and (roughly) ~ 7 1 . 
Generally, filled capsids have higher spring constants. Data for the filled capsids 
is for p = 0.06hie/R (internal pressure of phage 029). 




Fig. 4. Global response of elastic spherical capsids to local deformations. (A) 
For small indentations only the top of the shell is flattened giving rise to a linear 
force-distance relation. (B) At a critical indentation £ c ~ h/R the capsid partially 
inverts its shape. 



where AE is the change in total energy caused by the increase A£ 
in X. Both hexamers and pentamers were displaced. The indented 
hexamers discussed here all lie in the center of the facet spanned by 
their three neighboring pentamers. 

As can bee seen from Fig. |2]the elastic response of the capsid to 
the deformation strongly depends on the ratio of bending and elastic 
energy characterized by the FvK number 7. For small 7 capsids have a 
nearly spherical shape [12] and behave like homogeneous continuous 
shells. For 10 < 7 < 500 the (dimensionless) linear spring constants 
K\ of both hexameric and pentameric regions follow a square root 
law for sufficiently small £ < £0(7) 



F = K£ = 4V7£ , 



[4] 



see Fig. [3] The transition to the nonlinear regime (at £ = £ c ) takes 
place at smaller £ c the larger 7, see Fig. |2]\ and B. For small FvK 
number 7 ~ (R/h) 2 corrections to thin shell theory of order h/R 
become relevant. 

For 7 > 500 the capsid has a strongly faceted shape since the 12 
disclinations buckle. Due to this structural inhomogeneity the elastic 
response of pentamers and hexamers is different, i.e. pentamers gen- 
erally become stiffer than hexamers with increasing 7. Furthermore, 
the abrupt shape transformation shown in Fig. [2f only occurs for 
loading on pentamers. 

As Fig. [3] shows, the difference between the spring constants for 
pentamers and hexamers increases with increasing 7. At 7 « 2000 
the spring constants are already well separated explaining the exper- 
imentally measured bimodality of spring constants of 029 [9 |. From 
our numerical data we find K\ ~ 7 1 / 3 for hexamers and K\ ~ 7 1 
for pentamers. When loading on pentamers is parallel to ridges then 
Ki 7 s / 6 in agreement with the findings of Refs. 12711281 . 

The elastic behavior of spherical capsids under axial point force 
can be understood by simple scaling arguments l29l of an elastic 
sphere with energy E — E e + Eh. For small indentations X the top 
of the sphere is flattened in a circular region of diameter d, see Fig. 
|4]\. In the flattened region the radius of curvature ~ d 2 /X and merid- 
ians are compressed by an amount ~ £. In equilibrium one then finds 
d ~ R/j 1 / 4 yielding a total energy E/Kb ~ 7^ 2 £ 2 of the deformed 
sphere. Thus, the force-distance relation is linear with Ki ~ 7 



1/2 



In fact, a more rigorous analysis of the linear regime yields BOl 
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Fig. 5. Numerical indentation experiment on </>29 (with 7 = 1778). At an 
indentation of 0.13R the mechanical response of the pentamers suddenly de- 
creases due to shape inversion, see Fig. [4}3. Linear spring constants are deter- 
mined by fitting the region of small £, implying F = 258£ for pentamers and 
F = 139£ for hexamers. The inset displays the same curve as it would be mea- 
sured in an experiment with an elastic cantilever with spring constant K c = 100. 
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Eq. © for 7 > 10, see Fig. 1 in the supplementary material. For 
7 < 10 the scaling picture breaks down since the dominant bending 
energy is no longer concentrated in a small region around the pole. 

At a critical indentation £ c ~ /i/i?^7 _1//2 the shell undergoes 
a larger shape transformation in which a circular region in the upper 
part inverts its shape, see Fig. [4j3. Here, the deformation energy is 
concentrated in a ring of radius r and width d. The radius of curva- 
ture of the ring is ~ Rd/r and meridians are shortened by a factor 
~ rd/R 2 . By using again d to equilibrate bending and elastic energy 
(i.e. d ~ R/j 1 / 4 ) one finds E/ Kb ~ £ 3 ^ 2 7^ 4 leading to a nonlinear 
force-distance relation F ~ ^}^ 2 . 

For sufficiently large 7 the inversion transition of pentamers 
(shown in Fig. |4j3) is a first order transition (in £) as can be seen 
from the F(£) -curve in Fig. [2£. There, only the green lines corre- 
spond to conformations of minimal energy. The blue lines represent 
intermediate (metastable or unstable) conformations connecting the 
weakly and strongly deformed capsid shapes. For smaller 7 the tran- 
sition is continuous, see Fig. 03. 

With our numerical simulations it is also possible to extract pre- 
cise material parameters of experimentally investigated viruses. To 
do so we have simulated the conditions corresponding to the SFM 
experiments on 029 | 9 |. By using the FvK number 7 as fit parame- 
ter for the measured bimodality ratio Kf ent /Kf ex ~ 1.88 one finds 
7 c± 1778, in good agreement with the value estimated above from 
the geometrical parameters of Ref. (26). 

Then, Kb/R 2 — 1.14mN/m can be extracted directly from the 
dimensionless spring constant Ki = kiR 2 1 Kb shown in Fig. [3]by us- 
ing, e.g., the softest experimentally measured value of ki = 0.16N/m 
(with corresponding kf ent = 0.296N/m) |9|. The force scale is then 
Kb/R — 24pN. For empty CCMV the bimodality cannot be resolved 
experimentally (ki = 0.15 ± O.OlN/m ifTOl ), while numerically we 
find kf ex = 0.17N/m and kf ent = 0.12N/m. In this case, we extract 
K b /R 2 ~ 2.4 ± 0.2mN/m and K b /R~29± 2pN from the numerical 
data. 

With these material parameters a direct comparison with SFM 
experiments can be achieved. Fig. \5\ shows as an example the 
force-distance curves for hexamers and pentamers of 029. Within 
our numerical analysis it is even possible to take into account the 
elastic spring constant of the cantilever K c (measured in units of 
Kb/R 2 ). Then, a measured force F corresponds to a displacement 
X' /R — F/K c -\-^of the shell. Thus, as Fig. [5]shows for sufficiently 
stiff cantilevers the shape inversion of the viral shell leads to a dis- 
continuous force-displacement curve. The numerical results (shown 
in the inset) are in excellent agreement with the experimentally mea- 
sured force-distance curves |9 |. In particular, for both numerical and 
experimental deformations the shape instability occurs at £ c ~ 0.13 

Rupture. We have also analyzed the local response of elastic 
shells to loading by investigating the conditions under which rupture 
occurs. This process takes place in regions of high in-plane stress 
where bonds break due to over- stretching or compression. Thus, 
this analysis provides information about the interaction between cap- 
somers. 

As experimentally observed, phage 029 breaks at a critical in- 
dentation X r ~ 7.7nm ~ 0.35R |9|. To determine the associated 
change in bond length we have calculated numerically the shape of 
a shell with 7 = 1778 at this indentation yielding a maximal com- 
pression of c± 4.5%. The corresponding numerical rupture force is 
f r — 43Kb /R ~ l.OnN in agreement with experimental findings (9). 
Using this criterion, the numerics for empty CCMV predicts a rupture 
force f r c± 22Kb /R — 0.57nN while the experimentally determined 
value is f r ~ 0.63nN fTUl . 




Fig. 6. (Dimensionless) rupture force of an icosahedral capsid with (green 
curve) and without (red curve) volume constraint. Rupture is assumed to oc- 
cur at a critical bond length change of 4.5%. The rupture force increases with 
7 since it becomes harder to stretch/compress the material. The rupture force 
shows a 7" dependence. We find a = 0.65 ± 0.03 ~ 2/3 for low 7 and 
a = 0.53 ± 0.04 ~ 1/2 for large 7. The crossover takes places at the inver- 
sion transition shown in Fig. [4] The insets show ruptured structures for 7 = 10 
and 7 = 3162. The bonds where rupture occurs first are shown as red dia- 
monds. 

In the following, we analyze rupture of viral shells as a function of 
7. Fig. [6]shows the 7-dependence of f r . For 7 < 100 rupture occurs 
in the linear regime, see Fig. |4]\. Here, rupture is caused by compres- 
sion of meridians, see left inset in Fig. [6] In the linear regime, we 
observe that the rupture force scales approximately ~ 7 2//3 . However, 
most viruses have higher FvK numbers and will therefore rupture in 
the nonlinear regime after shape inversion, see Fig. |4}3. Here, rupture 
is caused by circumferential compression in the highly bent rim. In 
this regime, the rupture force scales (roughly) ~ 7 1 / 2 . 

We have also investigated the influence of the capsid geometry 
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Fig. 7. Rupture force in units of Kb/R plotted across the surface of spheri- 
cal and sphero-cylindrical capsids with 7 = 100 and 7 = 1000. Rupture is 
assumed to occur at a critical bond length change of 4.5%. 
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on rupture by comparing an icosahedral and a sphero-cylindrical cap- 
sid. This is motivated by the fact that even-numbered T-phages and 
029 have additional capsomers in the equatorial region elongating 
the icosahedral shape. For example, the head of T4 is an icosahedron 
with triangulation number 13 and one additional ring of hexamers 
||3T1 . The geometry of the sphero-cy Under used in our simulations 
is that of 029, i.e. an icosadeltahedron with triangulation number 3, 
elongated by two additional equatorial rings of radius R [26|. The 
force was applied perpendicularly to the axis of symmetry. In order 
to avoid tilting of the structure, the lowest 10% of vertices of the shell 
were kept fixed in the simulations. 

Fig. |7] shows the rupture force mapped across the surface for 
7 = 100 and 7 = 1000. In the corresponding numerical simula- 
tions every vertex was indented in steps of R/250 followed by the 
calculation of the new equilibrium conformation. When stretching 
of some bond exceeded 4.5% the applied force was determined. For 
7 < 7b below the buckling threshold the shells have a uniform elas- 
tic behavior across the surface. Since the force exerted on the caps 
has a tangential component, part of the displacement directly leads to 
bond stretching. Therefore, the caps are the regions of highest insta- 
bility. Above the buckling threshold, the pentamers are rigid cones 
but the space between them becomes more flexible. When pushing a 
pentameric region the displacement is mainly transformed into com- 
pression of the ridges. Pentamers are thus the most unstable region. 
On average a larger (scale-free) force is required to break the capsid 
for larger FvK number 7. 

The spatial distribution of deformation energy can serve as a mea- 
sure for the tendency to rupture. To mimic an ensemble of indentation 
experiments, we carried out simulations in which each vertex of the 
shell was indented by a constant £ (£ = 0.1). For all vertices the defor- 
mation energy was recorded for all deformations. In order to find the 
spatial distribution of deformation energy in the ensemble, the elastic 
energy per vertex was averaged over the ensemble measurement. Fig. 
[8] shows a map of the average spatial distribution of elastic energy, 
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Fig. 8. Averaged distribution of elastic deformation energy for an ensemble 
measurement where every vertex was indented by an amount £ = 0.1. In the 
simulations, loading was perpendicular to the axis of symmetry. The ensemble 
average shows that the elastic energy is concentrated in a small region just below 
the caps. The capsid will most likely break in this region. 
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normalized by the total elastic energy put into the system during the 
measurement on the ensemble. 

In our simulations we assume a uniform distribution of binding 
energies between capsomers. In reality, binding energies will show 
a spatial variation and an ensemble of indentation experiments on 
capsids will yield different rupture probabilities than the numerical 
data presented in Fig. [8] However, with a direct comparison of the 
two approaches conclusions about the distribution of binding poten- 
tials can be drawn. 

Filled Capsids. Phages infect their host cells by penetrating the 
cell membrane and rapidly injecting the DNA into the cell plasma. 
Therefore, their DNA is stored under high pressure. The pressure p 
of DNA inside 029 is of the order of p ~6MPa HIE) or in units of 
elastic parameters p ~ 0.06n e /R- 

Other viruses, like CCMV, self-assemble inside the host cell and 
enclose the DNA. Correspondingly, their internal pressure is much 
lower than that of phages. Typically, p =latm= 0.0lK e /R is a good 
estimate for the internal pressure of a filled CCMV capsid. 

To investigate the influence of internal pressure on the elastic prop- 
erties we have numerically simulated filled capsids with the methods 
described above. However, here the additional constraint has to be 
taken into account that due to the presence of DNA inside the capsids 
the enclosed volume is fixed. In order to mimic a packed 029 phage, 
first the equilibrium conformation under pressure was determined by 
minimizing the total energy with an additional pressure contribution 
E p pV with p 0.06K e /R. 

Experiments on full CCMV were performed in Ref. H01 . A 
direct comparison of our numerical simulations with these exper- 
iments yields (for p = 0.0lK e /R and volume V =const.) the 
spring constants fc z hex ~ 0.19N/m, kf ent ~ 0.12N/m, and a rup- 
ture force f* uU ~ 0.83nN. These values are again in good agree- 
ment with the experimental findings (ki — 0.20 ± 0.02N/m and 
f r ul1 = 0.81 ± 0.04nN). There are no SFM-experiments on filled 
029 phages to which we could compare our numerical simulations. 

Generally, the volume constraint leads to increased linear spring 
constants, see Fig. [3] This can be understood in the scaling picture: In 
order to preserve constant cross-sectional area, the local compression 
at the point of loading must be compensated by an expansion of the 
equatorial area. Therefore the influence of the volume constraint be- 
comes more apparent for larger 7, i.e. for shells dominated by elastic 
energy. 

At high pressures the circumferential stress at the fivefold discli- 



0.12 



0.115 - 

0.11 - 

,— , 0.105 - 
DC 

> 0.1 - 

°- 0.095 - 
0.09 - 
0.085 - 




°'°^e0 1e1 1e2 1e3 1e4 1e5 1e6 

Y 

Fig. 9. Rupture pressure (i.e. pressure at which the relative bond stretching 
reaches 4.5%) for an icosahedral shell as function of 7. 
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nations is balanced by the volume contribution pV. Thus, pentamers 
do not form rigid buckles and hexamers remain flat. For exam- 
ple, the DNA pressure of 029-phages reduces its aspherity a from 
a c± 1.4 • 10 -3 for an empty capsid to a ~ 0.6 • 10 -3 . This similar- 
ity between hexamers and pentamers is reflected by the fact, that the 
elastic response of hexamers and pentamers becomes similar for high 
internal pressure, see Fig. [3] 

The rupturing behavior of filled capsids is shown in Fig. [6] The 
rupture force is larger for filled capsids than for empty ones. Due to 
the internal pressure the non-indented shell is already stretched. This 
compensates the force-induced compression in meridional direction 
at the poles thus reducing their tendency to rupture. Rupture of filled 
capsids is caused by circumferential expansion at the equator. 

Osmotic Shock. One way to extract DNA from viral capsids 
is to put them under osmotic shock. Under these conditions some 
viral capsids (e.g. T-even phages) rupture [32 | while others (e.g. T- 
odd phages) stay intact. Thus, these experiments also elucidate details 
about capsomer-capsomer interactions and we have analyzed whether 
rupture induced by internal pressure and by external force are related. 

To do so we have numerically determined the shape of the cap- 
sid under internal pressure p and measured the bond length between 
neighboring vertices. It was assumed that at rupture pressure the criti- 
cal bond stretching is 4.5% as for rupturing due to bond compression. 
Fig. |9]shows the corresponding rupture pressure p r (in units of K e /R) 
as function of 7 for an icosahedral shell. 

As can be seen from Fig. |9] capsids with a buckled ground state 
(7 > lb) rupture at lower p r than unbuckled ones. One should note, 
that for rupture induced by an external force, rupture forces are gen- 
erally higher for 7 > 7^ than for 7 < 75, see Fig. [6] However, 
as can be seen from the inset in Fig. [6] in the buckled configuration 
the force-induced deformation propagates over a larger area than in 
the unbuckled state. Therefore, it is possible that the external rupture 
pressure for force-induced rupture is also lower for buckled capsids. 
But this questions needs to be addressed in a more detailed analysis. 

For small 7 the rupture pressure reaches a constant value, see Fig. 
|9] In this limit deviations from spherical shape are small and the in- 
ternal pressure simply leads to a form-invariant up-scaling of the viral 
shell. Thus, changes in bending energy can be neglected. By equi- 



librating the expanding pressure force and the restoring force arising 
from bond stretching for a triangulated sphere the relative change of 
bond length is found to be pR/ (3« e ). The numerically found value is 
somewhat lower caused by the small deviations from spherical shape. 
The same calculation shows that a triangulated sphere with a har- 
monic stretching energy is only stable for pressures p < 0.75ft e /i?, 
independently of 7, see supporting information. 

Summary and Conclusions. Recent SFM-measurements have 
revealed the elastic properties and mechanical limits of viral capsids. 
In this paper, we have shown that numerical simulations are a pow- 
erful tool in complementing these experimental efforts. We have ex- 
amined the reversible and irreversible mechanical behavior on local 
and global scales. In particular, we have shown that the elastic param- 
eters characterizing the mechanical properties of phages and viruses 
can be determined very precisely by a direct comparison of numerical 
and experimental data. With the numerical methods presented here, 
one is able to resolve the mechanical response of viral shells to an 
external force in high detail. In particular, since we are able to dis- 
tinguish the response of hexamers and pentamers we can identify the 
origin of the experimentally observed bimodality of the elastic spring 
constants. We also make predictions for the elastic response and rup- 
turing behavior of both empty and filled capsids as function of the FvK 
number 7. Furthermore, in our simulations the mechanical limits of 
viral shells can be probed in an ensemble where every capsomer is 
indented. The comparison of the corresponding rupture map of the 
shells with experimental data on SFM-induced rupturing offers new 
methods in experimentally probing the local protein-protein interac- 
tions. We can also use the rupture criterion to predict the maximal 
sustainable internal pressure of capsids. 

So far, we have focused on shells with the shape of an icosadelta- 
hedron. However, many viruses such as phage HK97 [ 33 1 or the 
animal virus polyoma [34] are chiral. The influence of chirality on 
the elastic properties and mechanical stability will be addressed in a 
forthcoming publication 1501 . 
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